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INTRODUCTION 

This document describes the results of a small project to 
extend and apply the analytical techniques (described in Ref. 1) 
which were developed under a previous NASA contract to predict 
the ignition and performance transients of segmented, solid pro- 
pellant boosters. In particular, the project was directed at 
developing techniques for analyzing the data from the first stc* c 
test of the Solid Rocket Booster (SRB) for the Space Shuttle (i.e., 
DM1) and predicting the performance of the second static test, 

DM2. 

The previous study (Ref. 1) developed analytical techniques 
for the type of segmented motor configuration and ignition events 
illustrated in Fig. 1. This was accomplished by accounting for 
(1) the temporal and spatial development of the flow field in 
segmented motors set up by the igniter discharge, (2) ignition 
and flame spreading coupled to the main chamber and slot flows, 

(3) the large velocity, pressure, and temperature gradients that 
occur during the early phases of ignition, and (4) the interactions 
that combine to produce peak pressures, (e.g., gas flow into and 
out of the slots, compression of chamber gases during pressurization, 
erosive burning, and mass-added effect of igniter discharge). 

The work to be summarized in this report has been catagorized 
as follows: 

1) Analysis and prediction of thrust dynamics during 
flame spreading and chamber pressurization. 

2) Extensions to account for the close coupling between 
the igniter and the head-end segment of the Space 
Shuttle SRB. 

3) General improvements and extensions to the computer 
program to facilitate parametric studies. 

As the pursuit of higher performance rocket motors continues 
and the present configurations arc upgraded (i.e., by increasing 
loading density, using high-performance propellant, and extending 
ambient temperature range) , the methods extended during the project 
will enable analysts to predict items such as ignition delays, rate 
of thrust increase and maximum pressu es with much greater confi- 
dence than was heretofore possible. 
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PYNAMIC THRUST DURING IGNITION TRANSIENTS 
Physical Situation Being Analyzed 

Rocket motor thrust is the net force exerted by the gases 
flowing in the rocket motor on the internal solid surfaces. Thrust 
can be calculated by either summing the resultant forces on the 
solid surfaces or by finding the change of momentum of the com- 
bustion gases in the chamber. 

In this discussion, a distinction is made between static 
thrust and dynamic thrust. Static thrust is the widely used con- 
ventional concept of thrust produced by gases flowing through a 
rocket motor nozzle. Static thrust is usually expressed as 

^noz - ^FAm^E, stag^t ^ 

where A. is the threat area, p_ . the stagnation pressure 

t f S'CclC^ 

at the nozzle end of the chamber, and c F ^ m is a coefficient 
which accou r .cs for the change of momentum as the gases flow from 
the chamber to the nozzle exit plane. The concepts of calculating 
C FAm so as to account f° r divergence losses, flow losses, flow 
separation, expansion ratio, etc., are well documented in books 

on propulsion. Static thrust is associated with quasi-steady 
flow; even though chamber pressure or throat area may be changing 
rather rapidly, the flow fields in the chamber and nozzle are usu- 
ally considered to be a series of steady states. Thus, static 
thrust is based usually on a series of instantaneous values of 
P E s t a g and A t* The conce P ts of dynamic thrust are rarely 
considered and not included in any of the standard texts or 
monographs on rocketry. 

Dynamic thrust is that component of total thrust which is 
produced by such events as rapidly accelerating changes in the 
internal rocket motor flow field, rapid relative movement of 
the propellant or inert hardware, sudden changes in flow area 
(e.g., nozzle removal or vent port opening), and ejection of 
objects through the nozzle. Dynamic thrust events noticeable on 
thrust versus time record are often referred to as thrust spikes. 
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Dynamic thrust is seldom of interest in the applications of 
conventional rocket boosters and sustainers. The only common 
rocket motor situation which produces sufficiently rapid changes 
in the momentum of the internal gases and discernible thrust 
spikes, is rapid ignition. In special cases, extinguishment by 
opening the case or by removing the nozzle can produce prominent 
thrust spikes. For most situations, the total impulses associated 
with these events are insignificant fractions of the total. 
Furthermore, because of limitations of thrust stand dynamic 
responses, thrust measurements during pressurization are usually 
very inaccurate (e.g., overdamped instrumentation to compensate 
for thrust stand ringing) -r.d, thus, thrust spikes during ignition 
are not observed. However, the boosters being developed for the 
Space Shuttle are sufficiently long that the dynamic thrust events 
associated with ignition and pressurization transients are dis- 
cernible. Indeed, the thrust versus time records of the first 
static te^t of the Space Shuttle Solid Rocket Booster (DM1) reveal 
clearly transient forces on the order 300,000 lbf during the first 
200 ms of motor operation. 

Figure 2 is a plot of the measured thrust versus time record 
from the first static test of the booster for the Space Shuttle. 

The several prominent events associated with the initial thrust 
are indicated. Coinciding with (what is calculated to be) the 
onset of flame spreading is a rapid increase in thrust. The 
thrust between 0.02 and 0.1 seconds occurs before flow through the 
nozzle is established and, thus, is attributed to the momentum 
changes that occur within the chamber. The periodic disturbances 
as thrust increases are attributed to wave action within the 
chamber and to oscillations of the propellant and case masses. 

The analysis summarized in this report is part of an effort 
to explain the observed thrust events. In particular, the analysis 
is directed at the effects produced by the developing flow field 
within the chamber. Another group is analyzing the interactions 
attributable to acceleration of the propellant and the inert parts. 
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Fig. 2 Dynamic effects noted on plot of measured net thrust 
versus time for the first static test of the booster 
for the Space Shuttle (DN1) . 



- 6 - 


General Form of Momentum Equation 


In this Section, rather general forms of the momentum rela- 
tionships are discussed for the purpose of defining those inter- 
actions which are pertinent. In the next section, a more manage- 
able set of relationships are applied. 

Let n (x,t) be any summable continuous function such as the 
density p, or momentum pV. The "content" of V is ndv. 

J 

The rate of change with time of that integral is given by Reynold's 
transport theorem. 


±_ f 

dt J 

V*(t) 


ndv = 


lt dv + 


v*(t) 


’ 

S*(t) 


nb 


ndA 


( 2 ) 


d/dt means: following the content within the bounding surface 

S*(t) which may move with £(t). 

9/3t means: fixing the references in space (relative to inertial 

frame of reference chosen) and changing with time. 

If the control volume chosen is a material volume [i.e., nc 
material element within S*(t) is ever crossing S*(t) and 
S*(t) is composed of the same material elements throughout the 
entire time of consideration] instead of d/dt in Eg. (2)/ write 
D/Dt to denote material volume and instead of b, u is used. 

So that 


D 

Dt 


ndv 

v*(t) 


3ru 

3t dV 


v*(t) 


f 

nu 

s*(t) 


ndA 


(3) 


The momentum balance equation for a material volume is: 


D 

Dt 


j 


pudv = 


V*(t) 

> V- ' V 

rate of change 
of momentum 


pSdv + 

V* (t) 

_ * 
body 
force 


$dA 

S* (t) 

1 11 - ■ 11 

surface 

force 


(4) 


Now, in the general case, a material volume is not used. Such is 
the case in thrust calculations. Thus, it is necessary to find 
the momentum equation in the most general cast when the control 
volume boundary moves with some b / u. 
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Using Reynold's theorems yields: 


5t f t,Sdu = 1 

v*(t) 


3t 


(pu)dv + 


(pu • u) • ndA 


V*(t) 


S* (t) 


(5a) 


d 

dt 


pudv = 


V* (t) 


g^(pu)dv + 


V*(t) 


(pu • b) • ndA 


(5b) 


S*(t) 


Though the control volumes' boundaries of Eqs. (5a) and (5b) move 
with different velocities u and b, at the time of investigation 
they momentarily coalesce, so that V*(t) and S*(t) in Eq. (5a) 
are identical with V*(t) and S*(t) in Eq. (5b). 

Summing Eq., (5a) and Eq. (5b) yields: 


D 

Dt 


pudv = ^ 


V* (t) 


V*(t) 


■+ f 
pudv + j 


pu 


(u - b) 


s* ( t) 


ndA 


( 6 ) 


Note: 


§t j p5dv * it 

V*(t) 


pudv / j dv 


V* (t) 


V*(t) 


The reason is that the volumes coincide only at time t, but not 
at time t + At, see illustration in Fig. 3. 

From Eq. (6) and (4) 


d_ 

dt 


pudv + 


• ndA = 


)($dv 


pu • (u - b) 

V*(t) S*(t) L VMt) S*(t) 

In the next section, Eqs. (6) and (7) are applied to calculate 
thrust transients during a static test. 


TdA (7) 


Motor in Thrust Stand 

As shown in Fig. 4, the rocket motor is free to move on its 
stand, exerting its thrust on the load cell which may move some 
distance e under the thrust reaction loading. The inertial 
frame of reference is naturally chosen as the load cell. 


t - time 

x - spatial location 



u and b are measured relative to some 
inertial frame of reference. 

(a) Illustration of parameters used in development of 
equations. 




■+— x 


V(t) 




j pu(t+At)dv 
V(t+At) 

r - 

pu (t) dv 
V(t) 


V(t) S V (t+At) 


pu(t+At)dv 

t) 


(b) Movement of control volumes. 


Fig. 3 General control volumes fur development of 
momentum relationship. 
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am 


ex 


(a) Overall view of motor attached to load cell 




x 


(c) Section through motor revealing stresses in case and 
propellant . 


Fig. 4 Control volumes and nomenclature used in development of 
dynamic and static thrust relationships. 
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The control volume is chosen to be the entire rocket motor. 
Now look at Eq. (7) and identify the various parameters appearing 
there. Since only the x direction will be considered, the 
vectorial notation will be dropped. 

The fc lowing descriptions apply: 

b - The control volume boundary velocity relative to the 
load cell. In this case it is the velocity of the 
rocket motor outer surface. If the motor case is 
considered to be rigid, then b is simply the rate 
at which the motor covers the distance e. 

u - The velocity of any material inside the control volume 
(relative to the load cell' . Conveniently, the veloc- 
ity of the gas inside the luotor is denoted as u^. The 
velocities of other components inside will be denoted 
as : 


G 

T 


Up r - propellant segments 

u - motor case segments 
case 3 


u / u / b 
pr case 

is very small , b 


u and u,,,^ << u . 
pr case g 

Note that u . u , 

g pr 


u 


case' 


It should be emphasized that due to elasticity 

and also u and u << u_. Since 

<< u 

g 

b ere functions of x and t and are taken positive 
if in the +x direction and negative if in the -x 
direction. 

Gravity, is zero in the x direction. 

The outer pressure force resultant in the x direction. 


Now look at each of the terms in Eq. (7) . The symbols 


V 


p pr' 


) are used for denoting densities of the gas, 

propellant, and case segments inside the control volume. In the 
sections of the control volume where there are no solid propel- 
lant segments, p =0. Where there is no motor case, p = 0 

pr Case 

and so on (e.u., in the port area section p = p = 0. Or. 

J case pr 

the outer shell, p = p =0.) dA in the integrals is taken 

A ^ A 

positive, if n is in the x direction, and taken negative if n 
is in the -x direction except in the surface force integral 
where it is taken always positive. Accordingly, the first term in 
Eq. (7) becomes 
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d_ 

dt 


I pudv = ~r— I ( p u + p u +p u )dv 
J M dt j M g g pr pr case case 


( 8 ) 


Since u >> b,u ,u , u can be taken simply as the value 
g pr case' g c * 

calculated from the set of flow equations inside the motor. Since 
ppr»Pcase >:> p g' eac ^ integrands in Eq. (8) may have signif- 

icant contributions and their relative magnitude should be evalu- 


dv 


ated before neglecting any of them. Generally, d/dt J p^u^ 
is the dominant term. 

Similarly, the second term in Eq. (7) becomes 

pu(u-b)dA = [p u (u -b) + p u (u -b) + p u (u -b) ]dA 

J * J lK g g g K pr pr pr K case case case 

S S 

(9) 

Obviously, the manner in which the control volume is constructed 
results in 


p u (u - b)dA = 0 
pr pr v pr 


(10) 


(i.e., no propellant on the boundary Pp r = 0) . 


p u (u - b) dA = 0 
case case case 


(ID 


(i.e., boundary is composed of motor case). 


p u (u - 

J g g g 


b) dA = p u A 

’ex ex ex 


( 12 ) 


(i.e., only at the exit plane p^ f 0 and, generally, u g >> b. 
So that Eq. (9) reduces to 


| pu(u - 


b) dA = p u A 

ex ex ex 


(13) 


Furthermore, considering the third term in Eq. (7), 


pGdv = 0 


( 14 ) 


since gravity is zero in the x direction. 
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Finally, the remaining term in Eq. (7) can be simplified. 


f 


TdA 


-p A + p A 
ex ex am ex 


+ F 


(15) 


Thus, the idealized thrust equation which accounts for 
dynamic thrust is: 


= d_ 
dt 


(p u +p u +p u )dv 
g g pr pr case case 


p u A 
ex ex ex 


P A 
ex ex 


P A 
r am ex 


(16) 


The integral on the right hand side of Eq. (16) is the 
dynamic contribution to thrust F^ ; the remaining terms are 
the conventional nozzle thrust usually expressed as F 

noz 

Numerically the above integral should be computed as follows: 
the momentary volume V(t) is divided into N stations and the 
increment of time is At. Then, 


d_ ’ 
dt . 
V 


(P u 

g g 


p u + 
pr pr 


p u )dv = 
case case 


- N 

l 

i=l 


ou+pu +p u 
g g pr pr case case 


dv. 

i 1 t+' 

-* f 


At - 


r- N 


y 

i^rg g "pr pr 


p u + p u + p 

~ /«f /r ~ TM" r 


u 


case case | ^ i| t 


d v. 


At 


(17) 


Since the values of u and u are not known at this 

pr case 

time and since the primary purpose of this analysis is to account 
for the gas flow effects, the calculations of dynamic thrust in 
the ignition transients computer program are carried out using 
a simplified form of Eq. (11) . 


~N 


“ N 

7 (p u A ^) .Ax 

— 

y (p u A . ) .Ax 

g g port'i 

J- _1 

t+At 

g g port i 


F dyn At 

The thrust associated with the nozzle is calculated using the 
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conventional nozzle thrust prediction procedures for quasi-steady 
flow. Thus, the net thrust is 

F ~ F dyn + C FAm P E,stag A t 


Discussion of Considering a Partial Control Volume 

The motor is now cut at the plane x - x. The main difference 
in the calculation is the considerations of the tensions (T 

case 

and T pr ) which may amount to an important part of the total F 

and should not be neglected unless u , p , b , u , u 

g,x *x x pr,x' case.x 

are zero. Following the reasoning as before we obtain: 

Segment A: 


F = 


dt 


(pu+pu +p u )dv+p u A + 

g g pr pr case case M g,x g,x port,x 


pA . -pA -T -T 

x port,x r am xx case pr 


(18) 


Segment B: 


T + T + ^ 

case pr dt 


(pu + p u + p u ) dv - o u A 
g g pr pr case case g,x g,x port. 


B 


-pA . + p (A -A ) + p A + p u A =0 (19) 

*x port,x y am xx ex ^cx ex ^ex ex ex ' ' 

It is easily seen that the summation of Eq. (18) and Eq. (19) 
yields Eq. (16) . 

During ignition period of a long motor, cross section x - x 

and p = p and, therefore, 

X alu 


can be found 

where 

u = 

p — p * 

o 

ii 

a 

. Thus 

r ex r am 

ex 

T + T 

case pr 

+ at j 

(p u 
’pr pr 


Substituting in Eq. (18) yields: 

F = tt I (p u -i-p u -p u )dv + ir- I (p u +p u )dv 
dt J g g pr pr 'case case dt J pr pr case case 

V, 
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If V ft >> V B , i.e., x - x is close to the nozzle, we can 
calculate F considering only Segment A and using a reduced 
Eq. (18). 

F dt | ^ P g u g + p pr u pr + p case u case^ v 
V A 

Generally, however, the thrust reaction F cannot be calculated 

considering a partial control volume like Segment A because 

T and T are unknown, 

case pr 
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EXTENSION TO COMPUTER PROGRAM 

In the subsections that follow, the modifications to the 
program and the new input and output are described in modular 
form. The descriptions follow the formats used in Ref. 1. The 
new computer program is designated as HVTSEG2 . Table 1 is an 
input data set for an example case. 

In the extended version of the program the subroutine TIMEST 
was eliminated and its function is now performed by MAIN. The 
subroutine FCT was eliminated since it dealt with a particular 
slot flow option which was found to be unrealistic, i.e., large 
pressure differentials between the main chamber and slots do not 
occur in the SRB. A multi-point interpolation subroutine, ITERP1, 
was added as part of the procedure for calculating thrust differ- 
entials . 

Although numerous changes and improvements have been made to 
the Ref. 1 computer program, the input data decks from the Ref. 1 
version of the program can be used without change. 

The following changes to Section 4 of Ref. 1 should be 
noted : 

pg. 47 - 


pg. 48 
pg. 49 


Delete XG as an input. The program sets 
XG = XE - Ax/2. 

Add to description of TIGN: Corrected by program 

whenever TPI ^ TOREF. 

Delete DDHC as an input. The factor multiplying 
the heat transfer coefficient is now CHC(N). 

The first two lines should be changed to: 
or by a tab’ e 

r 0 * ex Pl°p ,T pi - T 0,ref> ^ 


(4-2) 
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Ignition of First Segment in Space Shuttle SRB 

The large surface area associated with the head-end segment 
of the Space Shuttle SRB produces a rather unusual ignition sequence 
Once the head-end pyrogen ignites the head-end segment, the mass 
flux from the head-end segment dominates the mass flux from the 
pyrogen by more than 7 to 1. Thus, the head-end segment is 
effectively the igniter for the remaining segments. Accordingly, pre 
dictions of the start-up transients and parametric studies to 
tailor the start-up thrust versus time program, require that 
careful attention be given to the ignition of the head-end segment. 

In the current design of the Space Shuttle SRB, the nozzle 

of the igniter is very close to the star-point tips in the first 

2 

segment. It has been pointed out that treating the head-end 
region as a uniformly heated port (as in the Ref. 1 analysis) 
results in unrealistically long induction periods. A better simu- 
lation will be obtained by considering at least three zones of 
heating and ignition in the first segment: (1) the more intense 

heating, rapid ignition, and high rate of flame spreading of the 
star-point tips, (2) the increasing heating rate and flame spread- 
ing rate in the axial slots as the gas generated by the burning 
star-point tips augments the axial flow of hot gases, and (3) 
the conventional heating of the aft portion of the first segment. 
Thus, the net effect will be to reduce the predicted induction time 
and to extend the flame spreading interval for the first segment. 
Recall that the results of Ref. 1 showed that flame spreading 
rate (downstream of the first increment) is relatively insensi- 
tive to the propellant property changes, whereas time to first 
ignition is affected greatly by small variations in the ignition 
criterion and heat flux. Accordingly, an important extension to 
the computer program was to include the capability of treating 
more completely the specifics of the igniter discharge pattern 
and propellant surface geometry in the head-end segment. 
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7 November 1977 
HVTSEG 


Corrections for Changes in Initial Temperature* 

The program input should be prepared by considering the 
following inputs to be determined at the temperature, Tq re j: 

TIGN = T . 

ig 

TFREF = T , _ 

f ,ref 

RREF = r , 

ref 

TIGTAB & MIGTAB 
PDATA & RDATA 

When the input initial temperature (or ambient temperature) 

T pi differs from Tq re j/ the following corrections are made 

to account for the temperature differential, AT. = T . - T_ _ 

0 pi 0,ref 

W " T ig + (c pr /C ch )4T 0 

V T pi = T f,ref + (o pr /O ch )iT 0 

r@T . = r _ exp (a AT rt ) 
pi ref p 0 

TIGTAB@T pi = TIGTAB/exp(ir k AT 0 ) 

MIGTAB@T pi = MIGTAB exp(ir k AT 0 ) 

r@Tpi = f (PDATA, RDATA) exp (o p AT Q ) 

All of these corrections can be suppressed by inputting 
all of the values at the desired temperature, T pi' an< ^ 
inputting T 0ref = T pi . 

The burning rate corrections can be eliminated by 
inputting a p = 0.0. 

The corrections to the igniter mass flow rate table can 
be eliminated by inputting = 0.0. 


*It is assumed that the reader is familiar with nomenclature 
and input description of Ref. 1. 
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21 November 1977 
HVTSEG 


Special Input for Control of Flame Spreading 

Time of ignition and flame spreading rate after first igni- 
tion can be studied parametrically by invoking special inputs. 
Since the heat transfer coefficient at each station can be 
adjusted individually, the time of first ignition can be con- 
trolled. Also, the fraction of each axial station that is ignited 
initially can be prescribed. At each axial station*, during 
flame spreading 

b = + (1 ~ kfl^/^ff^ 

When flame spreading is complete, 

b = p 
*w 

The heat transfer coefficient at each station is the value cal- 
culated by the empirical equation multiplied by 

C h,n 


The following inputs can be prescribed for each axial station, 
from 0 to NDELX: 

CHC(N) C. A factor multiplying the heat 

' transfer coefficient for the 

purpose of providing a prescribed 
adjustment at each individual 
axial station. The input DDHC la 
constant multiplier for all 
stations) is replaced by C, . 

(1.0) ft,n 


BFKN) b fl/n 

DELTTF (N) At.. „ 

ff , n 


The fraction of the perimeter 
which burns following first 
ignition. 

( 1 . 0 ) 


Time required for flame to spread 
over remaining perimeter, i.e., 


the fraction 1 - b 
be zero . 

( 1 . 0 ) 


fl ,n* 


Cannot 


s 


♦Mote that in the extended version of the program the burning 
perimeter is printed out as zero until ignition occurs and then 
the time-dependent output of burning perimeter reveals how 
flame spreading within a single increment increases with time. 
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8 November 1977 
HVTSEG 


Special Inputs to Facilitate Parametric Studies 


The following inputs are used to extend the output of the 
program, simplify some of the inputs, and to improve the simula- 
tion of actual rocket motors. 


Computer Symbol 
Sumbol in Text 

GAMAN 

noz 

PISUBK u k 

TUN t. 

in 

DELTTF At p 

DFSDT At Ap 


Description 


Units 


Effective ratio of specific 
heats of all combustion 
products in the nozzle. 

(Y noz " Y) 

Temperature sensitivity of 
pressure and burning time of 
igniter, for corrections from 


0 ,ref 

( 0 . 0 ) 


to 


Pi 


Initial time at which solution 
begins. 

( 0 . 0 ) 


K 


-1 


s 


Thrust differentials between the 
times t and t - Atp are cal- 
culated at each printout interval 
for the purpose of examining 
thrust imbalance between two 
motors. 

( 0 . 0 ) 


s 


Thrust increases between the pre- s 
ceding time interval At^p are 
calculated at each printout inter- 
val for the purpose of examining 
the rate of thrust increase during 
pressurization . 

( 0 . 0 ) 
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7 November 1977 
HVT 


Special Input for Nozzle Closure 

The program can approximate the effects produced by using 
a nozzle closure during the chamber filling phase of ignition. 

The solution operates on the assumption that the nozzle closure 
does not completely block the flow at the throat. Thus, there 
mut,t always be a small vent area in the closure, e.g., 10% of 
the throat area. For practically all situations, a nozzle 
closure which blocks 90% of the throat produces pressurization 
rates which are comparable to a completely blocked throat. 

Since nozzle closures require a finite time to be carried away, 
the solution causes the throat area to increase linearly from 
the small vent area associated with the closure to the final 
throat area. 

The inputs for this option are : 

Gunits 

Computer Symbol CGS(modified) 

Symbol in Text Description 

POPEN When POPEN < 7776 . 0 , POPEN is atm 

considered to be the nozzle-end 
stagnation pressure at which the 
nozzle closure begins to carry away. 

(7777.0) 

NBL (1) Vent area in closure is At/NBL 

and (2) the number of At's required 
to open the discharge area from A t / 

NBL to A^, after the nozzle closure 
begins to carry away* ir NBL*BLANK4 . 

( 10 ) 

BLANK4 The number of time stept> required 

to remove tne nozzle closure is 
NBL*BLANK4 . 

( 1 . 0 ) 
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1 December 1977 
HVTSEG 


Special Summary Output 


The following time-dependent variables are output in 
tabular form at intervals NPRINT. British units are used 


in the summary. 



Computer 

Symbol 



Symbo l 

in Text 

Deeoription 

Unit 8 

TI 

t 

Time 

s 

PHPSI 

Pi 

Static pressure at head end of 
motor . 

,>sia 

PEPSI 

P E 

Static pressure at nozzle end 
of propellant grain (or blast 
tube) . 

psia 

PESTAP 

Pe, stag 

Stagnation pressure at nozzle 
end of motor. 

psia 

M 

“e 

Mach number at nozzle end of 
propellant grain (or blast tube) . 


FLBF 

F 

Nozzle thrust based on^y on 
quasi-steady nozzle flow. 

lbf 

FDYNLB 


Thrust associated with transient 
chamber flow upstream of nozzle. 

lbf 

FIMBAL 

^imbal 

Thrust differential between the 
times t and t - Atf 
(Based on nozzle thrust.) 

lbf 

FRATE 

AF in 
At AF 

Thrust increase during the time 
interval At^p. (Based on nozzle 
thrust. ) 

lbf 

FMIGLB 

m. 

ig 

Mass flow rate from igniter. 

lbm/s 

FNET 


Sum of nozzle thrust and dynamic 
thrust, FLBF + FDYNLB 

lbf 

The 

Summary output is set up to be displayed in three forms 


1. at the end of each time-dependent out (i.e., part of the 
standard listing) . 

2. as a deck for input to plot programs (i.e.. Unit 8 in 
subroutine ANS is designated SYSOUT = B in the JCL) . 

3. As a separate listing (see Table 2) following the 
standard listing (i.e., Unit 9 in subroutine ANS is 
designated SYSOUT = A in JCL) . 

The JCL is indicated in Table 1. To suppress punching of card 
deck, change SYSOUT=B to SY.SOUT=DUMMY in JCL. 


-W* DATA TO BE EXECUTED FOLLOWS G0.SY5IN DC 

■ //5o*f?09foot dd syson?=Al 

-^*G&rFTn«F00l dd sysout=bJ JCL ror Summary table 

//GO. SYS IN DD * 1 

0H1 1 7f'£C77 PLANE SPREAD CONTRO AT HEAD END ,D 

68AME 


TABLE 1 

input for example 
case for HVTSEG2 

YNAMIC THRUST CHECKOUT 


TflAX*.002 
TPFINT*. 092 
DELTA T*. 0001 ft 
NDELX=24 
PAH=. 84717 
0NlT=-2. 

AT* 23 25. 81 
XP=2. 67 
XP* 1329. 972 
RUFSUR=0 .01 
FKPR=C. 0011 
^=r... ROTP =1.756 8 
: CPF*. 3 

SIGP*. 0009930 
- RREF* 1.0540 


ORIGINAL PAGE IS 
OF POOR QUALITY 


PFFF=68.08 


OREXP = 0.33 

EBC-0, 

DE='145.64 
CH = . 97844 
ALFAD=12.31 
NI'TAn=19 


NAPDVX = 25 
DDHC=C. 94 
TPI = 298 .8 *1 
TIGN=2685. 4 
TPR EF=3353. 


H=28. 18 
T?SCRI=850. 
GABA =1.1360 
GAB AM =1.1348 ] 
PISUBK=0. 0027 J 
SDATA =3 
DELTTE=0. 00906 
DFSDT =0.010 I 


new 


new 


NBL=10 

BF1= 0. 1,0. 1,0. 1,27*1.0 
CHC= 1.4,1.3,1.2,1.1,26*1.0 

DE LTFF= 0. 010, 0.007, 0.0 05, 27 *1. 

NPNPXT=2 

NPN FXT= 0 

6ERD 


0 . 1 . 


0.0200 

18.11 

0. 0280 

46 . 95 

0.0440 

406.05 

0. 05 2 C 

484.35 

G. 0640 

514.35 

0.0700 

516.69 

0.0880 

510.06 

0. 1800 

427.47 

0.2400 

392.42 

0.2880 

371.00 

0. 3440 

348.43 

0.3560 

339.86 

-Q.J80Q 

/1 3. 75 


Input to control 
flame spreading 


< 2 ^ 



0.4240 275.18 


0.4600 

253.38 
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0*5080 

232.73 

. ... 

-i.; .* ~r. 

Z*cl« 1 (continue 


0.5560 

219.50 


- .T“ -* ’ 




0.5600 

0.0 

- 

. •. 

- . 



2.67 

2538 i 44 

1154.45 

1154.45 

31. 22*25i - 

5W1. 

2.21 

57.97 

2543.62 

1156.49 

1156.49 




113.28 

2561.7 

1144.95 

1144.95 




168.58 

27Q6.H7 

278 .93- 

278.93 - 

- - . 

- 


223.89 

2612.1 

131.18 

781.18 




- 279. 19 

2645. 

182*21 

183.31 

329.52 877. 

20151. 

1.5 

^ 334.5 

2658.85 

182.75 

182.79 

-- -_2_ _ :_v. : 



389.8 

2705.32 

184.38 

184.38 




445. 1 

2774.53 

186.72 

186.72 




S~. 500.41 

2844.42 

189. '6 

189.08 

• . ~ : — - 

r- 


flP 555.71 

2915.19 

191.4 

»1.4 

.. >.• ’ ‘-v^_ . _‘\T" J 



" 611.02 

3029.79 

195.13 

195. 13 

649.52 2890. 

20001. 

1.5 

— 666.32 

2757.93 

186.16 

186.16 

• - - ---• - - - -- 

- - . - 


721.63 

2720.09 

184.88 

184.88 




776. 93 

2789.3 

187.22 

187.22 





^-032.23 285U« --. - - r- 

1053.45 3224.03 201.28 201.28 

1108.76 3731.7 216.55 216.55 

■ ^ 1164.^ -4318.43 > 332 *9^ ^-- 232*9 - - -^~r--=*s=- =■■-- 

i :~3»fi^---;-3B*iai 7293.2 15349. 7?W4^ 9.1 

• : 1274.61 5404.39 Zifefe - '' — ^:;r - ^ - 

-•• 1 329. 972 - -5594 .-67 285.-15 -- 2*5 .15 - — - • • 

0.6808 0. 1002 
4.7651 0.4384 

68.08 1.054 .-. .- ■ 

= tun ■■ ■ - ■ ••• : '. 

' '. T8AX-.01 ' ■/ -". ^ ■• :, 

DELPAC-1. 3 ~ ' ••• - ' 

r, f,n n 

f,NA ” F ^flAX- 1 ORIGINAL PAGE IS 

TfUX«0.Q9 OF POOR QUALITY 

REND V :-. - ' A 

-■ RNAHF -■■■-'■ ------- - -'• - — -- 

TMAX-.12 

r>SLFAC=.5 

TPBIHT=0. 001 

SENT 
R N A 1 ” 

T i A X= . 3 ■ ■ - 

n T 'LFAr=2. 

6 END 
SNA.1 1 : 

TXAX=.4 
TPFIMT^.OOS 
DBIFAC- 1 . 

RDND 

RNA*!? 

TPHIMT=.01 
D8LPAC= 1. 

?f!AI*.5 

REND 
£ HA HE 

TH AX* 1 . 

- REMO 


ORIGINAL PAGE IS 
_2^_ OF POOR QUALITY 

Table 2 Special summary output table for example case for HVTSEG2. 

1 HVT ROCKET HuTOR IGNITION PREDICTION - S5GHBNT5D - NO? 77 

INCLUDES: SPATIAL 6 TISF DEVELOPHENT OP P,U,ST AND PlAflE SPREADING 

DBl 1 7DEC77 FLA HE SPREAD CONTRO AT HEAD END , DYNAHI C THRUST CHECKOUT 
**** SnHBARY OUTPUT 0? TIHF DEPENDENT RESULTS ***** 


TIP 5 

PHEAD 

PNOZ 

PC-NOZ 

HNOZ 

THRUST 

FDTN F 

IBBAL F 

RATE SIGN 

FNtT 

SEC 

PSIA 

PSIA 

PSIA 

- 

LBF 

LBF 

LBP 

LBP LB/C 

IFF 

0.009 2 

14. 

12. 

12. 

0.02 

0. 

0. 

C. 

0. 17.6 

• 

C.0020 

15. 

12. 

12. 

0.02 

0. 

74 12. 

0. 

0. 17.6 

7'» 

0.00«0 

14. 

12. 

12. 

0.02 

0. 

1636. 

0. 

0. 17.7 

1 6 1 n . 

0.0060 

14. 

12. 

i2. 

0.01 

0. 

1148. 

0. 

0. 17.7 

1148. 

0.0080 

14. 

12. 

12. 

0.01 

0. 

1260. 

ft. 

0. 17.8 

1269. 

0.0100 

14. 

12. 

12. 

0.01 

0. 

1464. 

0. 

0. 17.9 

1 A 

0.0120 

14. 

12. 

12. 

0.01 

0. 

2166. 

0. 

0. 17.9- 

2166. 

0.0140 

14. 

12. 

12. 

0.01 

0. 

2149. 

0. 

0. 18.0 

2149. 

0.0 16C 

14. 

12. 

12. 

0.01 

0. 

2163. 

0. 

0. 18. C 

2163. 

0.0180 

14. 

12* 

12. 

0.01 

----- 0. 

— 2202- 

0.- 

O— 18-1 

— . 2202. 

0.02CC 

14. 

12. 

12. 

O.OT 

^ • 0. 

2262. 

0. ... 

O. 18-3 

2262. 

C. 0220 

14* 

12. 

12. 

5,01 

9* 

2865. 

-V-:-- 6- - 

-i O. 25.6 

2865. 

0.0240 

14. 

- 12. 

12. 

0.01 

_ o._ 

3064. . 

_ . 

0. 32. B 

. 5064. 

C.C260 

14. 

12. 

12. 

0.01 

0. 

2864. 

0. 

0. 40.1 

7-64. 

0. 3280 

14. 

12. 

12. 

0.01 

0. 

3296. 

0. 

0. 48.6 

3296. 

0.0300 

16— 

12. 

12. 

0.01 

o.- 

7460.- 

— - 0.-- 

. 0. 93.7 

460. 

0.0320 

17. 

12. 

12. 

0.01 

, 0. 

12499. 

o.- 

• 0. 138.8 

12499. 

0.034 0 

18. 

12. 

12. 

0.01 

0. 

18309. 

o; 

0. 183.9 

18 30 9. 

0.0360 

20. 

12. 

12. 

0.01 

0. 

24627. , 

- _ - Q. - 

0,_ 2 29.0 

24627. 

0.G33 0 

21. 

12. 

12. 

0.01 

0. 

31141. 

0. 

0. 274.1 

3 141. 

0.0400 

22. 

12. 

12. 

0.01 

0. 

37730. 

0. 

0. 319.2 

377 . 

C.0421 

24. 

12. 

12. 

0.01 

0. 

44389. 

0 

— . 0. — 364.3 

44339. 

0. 044 0 

25. 

12. 

12. 

0.01 

0. 

51129. 

- • 0. 

0. 408. 1 

5112'). 

C .0460 

26. 

12. 

12. 

0.01 

0. 

554 53. 

0. 

0. 427.7 

55453. 

0. 0480 

27. 

12. 

12. 

0.01 

0. 

58766. 

0. 

_ O. 447.4 

58766. 

C. 050 0 

28. 

12. 

12. 

0.01 

0. 

61490. 

0. 

0. 467 1 

61490 . 

0.0520 

28. 

12. 

1 2. 

0.01 

0. 

63688. 

0. 

0. 485-8 

6368 8. 

0. 054 G 

3C. 

12. 

12. 

0.01 

0. 

63948. 

0. ... 

0. 490.9 

63943. 

0.0560 

30. 

12. 

12. 

0.01 

A 

V • 

6 3271. 

c. 

O* 495.9 

63271 , 

G.359C 

31. 

12. 

12. 

0.01 

0. 

62158. 

0. 

0. 500.9 

62158. 

C.C605 

31. 

12. 

12. 

0.01 

0. 

61156. 

0. 

O* 5 05.9 

61156. 

C. 0620 

34. 

12. 

12. 

0.02 

0. 

61193. 

0. 

0. 510.9 

611+ 1. 

0.0640 

37. 

12. 

12. 

0.02 

0. 

64839. 

0. 

0. 515.6 

64839. 

0.0660 

43. 

12. 

12. 

0.02 

0. 

72564. 

0. 

.0* 516.4 

72564 . 

0.0680 

56. 

12. 

12. 

C.02 

0. 

94464. 

0. 

' 0. 517.2 

94464. 

0.07 CO 

73. 

12. 

12. 

0.03 

0. 

132362. 

0. 

0. 517.9 

132362. 

0.0720 

92. 

12. 

13. 

0.04. 

0. 

185618- 

0. 

J>. 517.1 

185613. 

0.0740 

109. 

13. 

13. 

0.05 

0. 

241191. 

0. 

0. 516.4 

24119 1 . 

0.0760 

125. 

13. 

13. 

0.06 

0. 

237555. 

0. 

0. 515.6 

287955. 

C. 0780 

143. 

13. 

13. 

0.09 

0. 

341420. 

0. 

0. 514.9 

34142'. 

0.0800 

16 1. 

13. 

13. 

0.12 

0. 

397931. 

0. 

0. 514.2 

397931 . 

0.0820 

178. 

14. 

14. 

0.17 

5. 

449776. 

0. 

O* 513.4 

449776. 

0. 0840 

193. 

15. 

15. 

0.22 

0. 

487593. 

. . 0. .. 

0. 512.7 

487591. 

0.0860 

205. 

18. 

18. 

C.25 

0. 

510843. 

0. 

0. 511.9 

51.' 

0.0880 

216. 

22. 

23. 

0.26 

35699. 

521846. 

0. 

0. 511.1 

5575'i . 

0.090C 

226. 

3C. 

31. 

0.26 

57143. 

521546. 

0. 

0. 509.3 

57R6 

0.C910 

220. 

36. 

37. 

0.26 

73947. 

446518. 

0. 

508.4 

520464. 

0.0920 

232. 

45. 

07. 

0.26 

101068. 

407510. 

0. 

0. 507.5 

508573. 

0.0930 

235. 

60. 

62. 

0.26 

145021. 

357513. 

0. 

0. 506.6 

502534. 

0.0940 

23e. 

82. 

85. 

0.26 

216007. 

306944. 

0. 

0. 505.7 

527951. 

0.0950 

242. 

115. 

119. 

0.26 

3244 00. 

260417. 

0. 

0. 504.8 

53481 ■ 

C. 0960 

246. 

158. 

164. 

0.26 

472083. 

214675. 

0. 

0. 503.9 

68675 3. 
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Table 2 (continued) 


, -C.097Q 250. 205. 213. 8. 0. 503.0 

0.0980 250. 245. 25*. 6.26 785125- 59012. 0. 0. 502.1 

^ 0.0990 - 258. 269. 279^ 8vif6- ^ 7 9 7 8 7. —9 6 3 3 9. - -frv O. 501.2 

0.1000 262. 275. 285. 0.26 9C4994. -103972. 832361. 0. 500.1 

0.1010 265. 269. 279. 0.26 882090. -114840. 783147. 808160. 499.4 

; 0.1020 269* - 258* 268. Q .26 6376 12 r -92069* 696229* 736772 . 498.5 

.= 0. 1 030 ... 273 . 247. 256*: 6i26 ^EMl?2.^r 57287. 580590 . 696196. 497.6 

_ 0. 1040 276. 237. 246. 0.26 7548«. -20005. 438109. 538869 . 496.7 

^ 0. 1050 280 . 229 . 238; 8 *-26 727995. 13220. 2654 36 . 403695 . 495.8 


0. 1060 283 . 225 . 234 . 0.26 712264. 38921. 82323 . 240303 . 494.9 

0.1070 287. 223. 232. 0.26 706443. 53586. -71297. 67086. 494.0 


15=6- 1080 290* 224^ 232. 0.26 797099*- 61686.-168923, —77933. 493. 1 

iMO. 1090 294. 225. 233. 0,26 711502. 66470* -If 3333. -tt«39. 992.2 

0.1100 297. 227 . 236 . 0.26 720562. 676 14- -16 3526. -184431. 491.3 

- 8.1110 300 . 230 . 239 . 0.26 731817* 66316.-108755.-158305. 490.4 

0.1120 303. 233. 242. 0.26 743140. 65667. -50569. -94707. 489.4 

0.1130 307. 237. 246. 0.26 754321. 65808. -2419. -36883. 488.5 


1140 3 18. 240. 

MgtlSO 313. 243. 

ris^ii&o 315.-3*8. 
---8*1 170 -319. -249. 


249. 


6 

T 




0. 1180 

323. 

253. 

262. 

0.26 

616002. 

62166. 

104900. 108905. 484.0 

0. 1190 

326. 

256. 

266. 

0.26 

830C19. 

59855. 

110069. 118523. 483.1 

-i 6.1200- -329. 

266* 

276*= 

0*26 

- fti ft ? n? 

57313* 

113076* 123648. 482. 2 



Note: There are three thrust outputs 

• - FNET = THRU ST -■» -TO YN 



where THSDST is 
and FD YN Is dynamic thrust. 


nozzle throat 


. . ' — - ' , oRS Qg ttL PAGErl S : ^ - 

OF POOR QUALITY 


784 3*56 . 
83M n. 

813461 . 

I < • ' . 

761 '50. 
745 74 9. 
733195. 
7348C1. 
741215. 
751185. 
760'' *0. 
768785. 
777972. 
788176. 
798127. 
808807. 
820129. 
-831 700 . 
843234. 
8545 * 5 . 
866193. 
878168. 
889 Q 75 . 
901515. 
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RESULTS FROM EXTENDED PROGRAM 

Two sets of calculations were performed for the purpose of 
demonstrating the features of the extended program: 

Case 1 - DM1 17DEC77 FLAME SPREAD CONTROLLED 

Input approximates situation for DM1 and the new input 
parameters (for ignition and fj.ame spreading) were 
selected to accelerate initial ignition of the first 
segment and delay the flame spreading along the longi- 
tudinal slots of the first segment. Input values are 
listed in Table 1. 

Case 2 - DM1 6JAN78 FLAME SPREAD NOT CONTROLLED 

Input is the same as Case 1, except the controls over 
ignition and flame spreading as imposed by the new 
input parameters are not invoked, i.e., 

BF 1=30*1 . 0 ,CHC=30*1 . 0 ,DELTFF=30*1 . 0 

Figures 5 and 6 are normalized plots of the igniter mass flux, 
pressure rise and thrust evencs during the induction, flame spread- 
ing, and initial pressurization phases. In Case 1, ignition of 
the starpoint tips occurs more rapidly than the uniform ignition 
situation in Case 2; thus, head -end pressure begins to rise rapidly, 
about 0.008 s sooner, however, as shown in Fig. 7, the average 
flame spreading rate of Case 1 is about 30% slower than Case 2. 

This may not be a result of the imposed delay in the flame spread- 
ing rate of the longitudinal slots, but rather a result of the down- 
stream region not being preheated by the 'gniter mass flow. How- 
ever, the slower flame spreading rate of Case 1 does not diminish 
the strength of the axial pressure wave advancing down the port 
nor does it diminish the dynamic thrust contribution. 

As shown in Figs. 5 and 6, the dynamic thrust contribution is 
prominent during the time of flame spreading and decays rapidly as the 
longitudinal pressure wave encounters the aft-end of the motor. 

Indeed, a sufficiently strong longitudinal pressure wave can produce 
a negative F^ as the wave encounters the aft closure. Figure 8 
reveals that the time relationship between the igniter flow arriving 
at the nozzle and the surge of hot gases from the head-end segment 
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0.60 0. BO 
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Igniter mass flux, pressure rise, and thrust 
events for CASEl 17DEC77 FLAME SPREAD CONTROLLED. 
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CASE 1 FLAME SPREAD 

IN FIRST SEGMENT CONTROLLED 

CASE 2 FLAME SPREAD __ 

IN FIRST SEGMENT NOT CONTROLLED 
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Fig. 
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Flame spreading for CASES 1 and 2 showing that 
the slower flame spreading rate in the first 
segment significantly decreases the average flame 
spreading rate over the length of the motor. 



Pressure vs time at five axial stations 
showing how flame spreading rate in the 
first segment alters the flow field 
development. ' ' 


Fig. 8 
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affects the magnitude of the thrust overshoot. Apparantly, if the 
igniter gases lead (by a sufficient period) the gases from the 
head-end segment, the F^ contribution is diminished and the 
effect of the strong longitudinal pressure wave is attenuated. 

The dynamic thrust component has the proper trends with respect 

to the observed p and F vs t for DM1. However, net thrust measured 

for DM1 is the resultant of fluid flow dynamics as well as the 

dynamics of the motor case axial expansion and axial movement of 

the large propellant masses. Thus, the values of F . for Cases 1 

net 

and 2 cannot be compared directly to the measured thrust values. 
After the Thiokol/Wasatch personnel complete their calculations of 
the case and propellant dynamics, the net resultant forces at the 
load cell can be approximated. 

The important conclusions from the Case 1 and 2 calculations 

are: 

1) The time period of the dynamic thrust coincides with the 
initial thrust rise and decay measured for DM1. 

2) The overall flame spreading rate and dynamic thrust con- 
tributions are significantly affected by ignition and 
flame spreading of the first segment. 

3) Longitudinal wave action within the chamber can produce 
accelerations in thrust similar to those noted in Fig. 2. 

4) The modified program is an effective means of evaluating 
methods of tailoring the rate of thrust increase during 
the start up transients. 
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A 


(t) 

C FXm 

F 

fc(x,t) 

A 

n 


Nomenclature for Analysis of Dynamic Thrust 
area 

local boundary velocity which may vary over the surface S 
thrust coefficient with losses taken into account 
thrust 

body force (like gravity, electric field, etc.) 
outward unit normal vector to the surface 


P 


pressure, static pressure unless designated as 
stagnation pressure 


S(t) = bounding surface 


T = force 

T(n,x,t)= surface force. Tensor, depends on orientation of the 
surface element defined by n (generally pressure, 
friction or thrust reaction in a test stand) . 


t = time 


u(x,t) = velocity of any material element inside V (may be 
a gas element, a motor case element, or a propellant 
element) . 

V(t) = a moving volume (not necessarily a material volume) . 
v = control volume 

x = axial distance from head end 


e 

ri 

P 


= movement of head end of motor 
= any summable continuous function 
= density 
Subscripts 


am 


ambient 


case = case 
dyn 


ex 


dynamic 

exit plane of nozzle 
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g 

stag 

t 

P r 

port 

x 


gases in rocket chamber 

stagnation 

throat 

propellant 

port area 

axial direction 
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